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We investigate the fluctuation effect of the di-fermion field in the crossover from Bardeen-Cooper- 
Schrieffer (BCS) pairing to a Bose-Einstein condensate (BEC) in a relativistic superfluid. We 
work within the boson-fermion model obeying a global U(l) symmetry. To go beyond the mean 
field approximation we use Cornwall- Jackiw-Tomboulis (CJT) formalism to include higher order 
contributions. The quantum fluctuations of the pairing condensate is provided by bosons in non-zero 
modes, whose interaction with fermions gives the two-particle-irreducible (2PI) effective potential. 
It changes the crossover property in the BEC regime. With the fluctuations the superfluid phase 
transition becomes the first order in grand canonical ensemble. We calculate the condensate, the 
critical temperature T c and particle abundances as functions of crossover parameter the boson mass. 



INTRODUCTION 



Fermion pairings are mechanisms for superconductivity and superfluidity. In the case of weak attractive interaction 
between fermions the pairings are well described in Bardeen-Cooper-Schrieffer (BCS) theory where fermion pairs 
are typically of a size much larger than the mean interparticle distance. In some sense the pairing can be regarded 
as taking place in momentum space. As the attractive interaction gets strong enough fermion pairs become real 
Q_i' bosonic bound states. Below a critical temperature macroscopically large number of these molecular bosons occupy 
JL ! the ground state and form a Bose-Einstein condensate (BEC). The idea of covering BEC regime in an extended BCS 
theory was proposed by Eagles Q and Leggett 0| • A quantitative description of the crossover at finite temperatures 
from weak to strong coupling regime was given by Nozieres and Schmitt-Rink Q , who wrote down the universal pair 
'■ wave function which can be reduced to the correct ground states at BCS and BEC limit. The BCS-BEC crossover 
' have been extensively investigated in cold atom systems since such a crossover can be observed in experiments by 
tuning the magnetic field around Feschbach resonance to achieve situations with different scattering length among 
fermionic atoms [!, 0, 0, S, El ■ 

Recently there is a growing interest in extending the theory of the BCS-BEC crossover to relativistic systems. The 
crossover from pion condensation to Cooper pairing of quarks and antiquarks at large isospin densities provides such 
an example [Tol. [Til [ill. [T3| . Another example is color superconductivity where quarks form Cooper pairs on the Fermi 
surface due to an attractive interaction mediated by gluon exchange [3, [TH, [Tg, [13, EE OH 0, E], H3| (for reviews, 
see, for instance, [2J, |25j, [26|, [27|, \2Q, |29j, [SO, [31()- Because of asymptotic freedom, color superconductivity at very 
large densities can be studied in a weak-coupling approach (20l . ED . \2% . I23L I32L l33l . I34I l35T | using perturbative QCD 
(see, for example, [36. l37l l3§10 . For moderate densities more phenomenological models such as the Nambu-Jona- 
29l. l39l] or self-consistent Dyson-Schwinger equations [13, Ell are needed. All these approaches 
usually assume BCS-like picture. When the coupling strength turns even stronger, the BEC-like diquark pairing 
may set in. A lot of authors used the NJL model to study the BCS-BEC crossover in relativistic superconductors 



or superfluids [42l, |43j, |44j, |45j, |46j, |47|, |48j- In our previous work we set up a theory for the crossover in relativistic 
superfluids which includes both bosonic and fermionic degrees of freedom [49]. It is inspired by the boson-fermion 
model of superconductivity [Bol. l5ll [H3| . We only addressed the mean-field approximation of the model. The crossover 
is realized by tuning the difference between the boson mass and boson chemical potential. 

In this paper we extend our previous work [49] by systematic incorporating the fluctuation effect in the crossover 
from propagating modes of the di-fermion field within Cornwall- Jackiw-Tomboulis formalism (CJT) [53]. Our starting 
point is the boson-fermion model with a global U(l) symmetry. Similar efforts have been made in the NJL model 
[45l [T^l . The CJT formalism has been employed to analyze fluctuation effects in color superconductors [HI, H3|- The 
inclusion of fluctuation contributions leads to the change of the second order phase transition of normal to superfluid 
to a first order one. The paper is organized as follows. In Sec. II we re-write the boson-fermion model in terms 
of the Higgs and Nambu-Goldstone fields. In Sec. Ill the effective potential is evaluated up to two-loops including 
the mean field as well as the fluctuation contributions. The Dyson-Schwinger (DS) equations are derived from the 
effective potential in Sec. IV. The gap and charge density equations are obtained in Sec. V by taking derivatives of 
the effective potential with respect to the diquark condensate and the chemical potential respectively. The numerical 
results are presented and analyzed in Sec. VI. Finally we discuss about our results and draw some conclusions in the 
last section. 

Our convention for the metric tensor is g^ v — diag(l, — 1, — 1, — 1). Our units are H = c = k B = 1. Four- 
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vectors are denoted by capital letters, K = K^ 1 — (fcojk) with k = |k|. Fermionic Matsubara frequencies are 
u> n = iko = (2n + 1)ttT, while bosonic ones are w„ = iky = 2nirT with the temperature T and n an integer. Our 
convention for the Feynman rules in finite temperature and density field theory are chosen as compared to those at 
zero temperature and density, where the propagator and the vertex for a quantum field 4>{X) {X is a real time-space 
coordinate) in the theory are like, e. g. iG(X 1 ,X 2 ) = (T<j){X 1 )(f>(X 2 )) and iT = i8 n /[8ct){X l )5(j){X 2 ) ■ ■ ■ 5ct){X n )]. 
At finite temperature and density, we replace iG — > — G and iT — » Y. The energy unit is chosen to be the fermion 
momentum cutoff A appearing in fermion momentum integrals. 



II. BOSON-FERMION MODEL 



In order to describe the fluctuation effects from the propagating modes of di-fermion bosons, we start from the 
following Lagrangian which respects the global U(l) symmetry corresponding to total particle number conservation, 

= -]fS^ + \{d v -ifi b 5 M )<S>\ 2 -m 2 b \<5>\ 2 

+- ($ f Wr * - $lf t*) . (i) 

In the first term or the fermion sector, \& = ( "f ] and '5 = (ip,ip c ) are spinors for fermions in the Nambu-Gorkov 
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space the inverse fermion propagator Sq 1 is given by 



(NG) basis, where the charge conjugate spinor is defined by ipc = Cip and ip c = i/j t C with C = ij 2 j°. In momentum 
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^° (r > ~ [ 7^ - M7° - m J ' [Z) 

where fi and m are the chemical potential and the mass of fermions respectively. In the second and third terms of 
the Lagrangian JTJ) the boson chemical potential and the mass are denoted by fi b — 2/i and m b . In the boson-fermion 
interaction part, the last term of the Lagrangian (HJ, we have defined the interaction vertices 

f = 2«7i 75 (J [J), ^ = -2^75(0 j), (3) 

with g the Yukawa coupling constant. We can write 

$ = (p + (f Q = -^=ifR + i( Pi) + Vo, 

where ipo is the expectation value of the vacuum or the zero mode of the boson field and Lp R I denote the real and 
imaginary parts of non-zero modes. With tpo the U(l) symmetry is spontaneously broken. To make our approach 
transparent we did not include in Eq. HJ the self- interaction of <I> such as <fr 4 , which does not change the results 
qualitatively. The quartic or higher terms in A (or <po) are present automatically from the effective potential. The 
Lagrangian {1} then becomes 

£(A )¥> ,*) = -I^- 1 ^ + ^I_Zli|A| 2 + ~ ^6)^| 2 - |V^| 2 - TT^gl^l 2 + I (cp+Wf vl/ - ^^ft^ 

1 5?o-i,t, 1 / i b- m 6iAi2 !/ _ > n -if <Pn\ , ,^\mf^R 



The action is /(A, <p, = J x £(A, ip, VP). The boson-fermion vertices in the second equality are defined by 

f R = V2(f-f f ) =i2V2g l5 a? G , 

f/ = -iV2(f + f f ) = -t2V2g Tj(J ? G , 

where cr^p are Pauli matrices in the Nambu-Gorkov basis. The inverse fermion and boson propagators S" 1 and D^ 1 
are given in momentum space by 

s _! _ _( 7 M P^ + /i7° - m i-y 5 A 
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+ n\ - m 2 b 2p b ip 
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where A = 2gtpv is the condensate. We have dropped the mixing terms of zero and non-zero boson modes since they 
vanish when carrying out the path integral. We see in (@]) that the di-fermion boson is described by two real scalar 
field. 

The Lagrangian ([I]) can also be expressed in terms of the Higgs and the Nambu-Goldstone fields by parametrizing 
the complex boson field as <I> = + i]o)e 2ie with rj + rjo = y/2\$\ and 29 = arg($). Here the phase field 9 is the 

Nambu-Goldstone field and r\ the Higgs one. We can choose the unitary gauge by rewriting the fermion field in the 

form and ^ = _ ig , I . Inserting the above expression for <& and \& into the Lagrangian ([1]), we obtain 



C(A,<p,V) = -\^S-H + ^^A^ + \[{d vV f + ^l-ml)^ 

1,„„^ . 1 



2 4g 2 2 

+2(r, + m ) 2 [{d v 6f - HbdoO] - ^(^fl)^37/»* + ^*T*, (6) 

where S~ x is given by Eq. ((5j) but with A = \/2g?7o which is different from the case of Q. Note that in the Lagrangian 
§6§ we still use the same symbol \& and 9 but now without the phase 9: '5 = ^ ^ ^ and \& = (ip, ipc)- I n the second 



to last term 03 is the Pauli matrix in the NG basis. In the last term we have defined a new vertex 

f = V2«37 5 



1 

1 

We can see that the mass term of 9 is absent in the tree level as expected for the Nambu-Goldstone mode. 

A few remarks about the Lagrangian (J4j) and |(6]) are needed. We see that the complex boson field ip in Q is 
transformed into the Higgs and Nambu-Goldstone fields in the Lagrangian (J6J . Both §4§ and (J6J are parametrizations 
of the same Lagrangian (U) which are equivalent to each other. We choose Q in this paper as it is more simple and 
transparent than ([6]). For example, in (j4]), in addition to the 2PI diagrams from the fermion and Higgs fields, there are 
also diagrams from the fermion and Nambu-Goldstone fields from ([6]), which are much more complicated. Another ex- 
ample is that the only two-point Green function for bosons from Q is (Pip(X)^ (Y)) ~ (Pry(A)e 2ie(x) r;(F)e _2ie(y) ), 
where P denotes the path-ordered operator, contains the correlations of both the Higgs and Nambu-Goldstone fields. 

Note that we used the roman letter S to denote the bare fermion propagator but with the condensate A, which is 
also called the tree level propagator. We will use calligraphic letter S to denote the fully dressed fermion propagator. 

III. CJT FORMALISM 

We start from the Lagrangian |(4]). The CJT effective potential reads [H3|, 

T(A,V,S) = -1(A) + - jTrlnlT 1 + Tt(D~ 1 V- 1) 

-TrlniT 1 - Tr^S - 1) - 2T 2PI (V,S)} , (7) 

where T>, S are full boson and fermion propagators and I{A) = /(A, 0, 0) is the tree-level action with all fields replaced 
by their expectation values (we have assumed (ip) = (SP) = ('J) = 0). The factor 1/2 for the fermion term — TrlnS 
comes from the doubling of degrees of freedom in NG basis. Z? -1 and S -1 the inverse tree-level propagators for bosons 
and fermions, 

= (8) 

« ^(x)^-(r)' w 

with i,j = R,I and S~ x is given by Eq. ©. Note that all propagators differ their normal forms by a sign, i.e. 
S — > —S and D — > — D etc. The contribution of 2-particle irreducible vacuum diagrams T2 is truncated to the 2-loop 
level as shown in Fig. [H 



T 2PI (V,S) « --^[VxyTt^xSxyTySyx]} , (9) 



where X, Y denote indices of color, flavor, space, Nambu-Gorkov and/or bosonic species, and the trace is of functional 
sense running over the general indices. The full propagators are denoted by calligraphic letters. The negative sign 
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Figure 1: (a) The 2PI vacuum diagrams up to 2 loops, (b) The boson self-energy derived from the 2PI vacuum diagrams, (c) 
The fermion self-energy derived from the 2PI vacuum diagrams. The solid and dashed lines denote fermion and boson full 
propagators. X, Y label indices of the color, flavor, space and/or bosonic species. 
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comes from the quark loop. A factor of 1/2 is from the Taylor expansion of the interaction and the other one is from 
the NG basis. 

The Dyson-Schwinger (DS) equations for fermions and bosons can be derived by taking derivatives of the effective 
potential with respect to propagators, 
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which lead to the DS equations, 
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S- 1 -E(X>,S). 



(10) 



(11) 



where the self-energies for bosons and fermions II and £ come from T2P1 in Eq. © and are shown in Fig. QJb) and 
(c) respectively. They are expressed by, 



iTr[fsfS], 



n(s,x>) 

E(X>,S) = Tr[OT5r]. 
Substituting the DS equations ifTTj) into Eq. we arrive at 

T(A,V,S) = -I{A) + ^{Tr\nV- 1 + r Tr[D- 1 V-n{V,S)V-l] 

-Trln,?- 1 - Tr[5 -1 5 + £(£>,S)S - 1] - 2T 2 pi{V 1 S)} 
= -/(A) + - {Trln [D-\l + DIP)] - Trln [5 _1 (1 - ST,)] - Tr(S5)} 

w -/(A) + - {TrlnZ?" 1 - TrlnS -1 + Tr[DU{V,S)}} . 



(12) 



(13) 



Hereafter we use the thermodynamical potential fi to replace the effective potential r, ft = Oq— T2PI, where the sum 
of first three terms in Eq. (fT3|) are denoted as Q,q, Note that the expansion in terms of self-energies made in the last 
line somehow leads to a partial loss of self-consistency, but it much simplifies the numerical calculation. One can also 
make an expansion in terms of the condensate A near the critical temperature where A is small [H, HH]. 



IV. PSEUDOGAP AND RENORMALIZED BOSON MASS 



In this section we address the DS equations ifTTj) for bosons and fermions. First we address the DS equation for 
fermions. When obtaining S from elements of S^ 1 , we will encounter a term (iS _1 )i2[(iS _1 )22] _1 ('5 _1 )2i which is 
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Figure 2: Diagrams for the gap and pseudogap. The numbers 1 and 2 are the NG indices. 
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illustrated in Fig. [2|a) . We will show that the correction to the fermion self-energy due to fluctuations of di-fermion 
bosons in Fig. [2](b) contributes to the diagonal parts of the full inverse propagator S^ 1 . 
The contribution from Fig. ifJb) can be written as 

En(X>,S) « {^[Aj-fiSf,-]} 

d 3 q 1 . _„ 1 



7oA"f ? — — „ -= (14) 



(2tt) 3 (po + Mb ) 2 - £ p 2 ,u p+ Vo + 9o - fi - eE p 



where we have taken the bare propagator S and D in |(5]) as approximations to the full ones. Here i,j — R,I denote 
the real and imaginary components, go is the bosonic Matsubara frequency go = i2nnT. We also uses the equalities 
Drr(Q) = Du{Q) and Dm(Q) — —Dri(Q), see App. ([A]). The energy projector is defined by 

A e = - (\ + e - 01 ' P + 7 ° m 



v 



where e = ±. The presence of En/22 is equivalent to (5 1 )n/22 gets a correction, where S 1 is given by Eq. J5|). We 
can verify it by seeing 

S n = [(5- 1 ) 11 -(5- 1 ) 1 2((5- 1 ) 2 2)- 1 (5- 1 )2i] _1 

= [(On - E n - (5- 1 ) 12 ((5- 1 ) 22 - E 22 )- 1 (5- 1 )2i] _1 

« [(^- 1 )ii-E 11 -(5- 1 ) 12 ((5- 1 )22)- 1 (5- 1 )2i] _1 (15) 



where En is assumed the value in Eq. lfl4|) . Note that in reaching the last line of the above equation we have 
neglected E 2 2 for clarity since it is sandwiched by (S' _1 )i2 and (5 l_1 )2i and complicates the calculation. We do not 
include fermion self-energy contribution to En, ~ Toff 2 &o ha %r [H, [HI coming from long range gauge interaction 
such as magnetic photon or gluon exchanges, since the gauge field is absent in the current version of our model. We 
can easily evaluate 

d 3 n . „ 1 1 



E 11 -(5- 1 )i2((^ 1 )22)- 1 (5- 1 ) 2 i « -16.g 2 T^ J t^ToA 



(2tt) 3 ,u p po + go - It - eE p (p + fi b ) 2 - B% 

where we have assumed that g <C p so |p + q| ~ p. The so-called pseudogap A pg can be defined [H, H3, H§| , 
A 2 1 , 2 T y f d 3 q 1 _ 2 f d 3 q 1 + f B (E b q - fi b ) + f B (E b q + fi b ) 
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d s q f B (E b q -p b ) + f B (E b q + M6 ) 
165 / (2^ 2E\ ' (17) 



where in the last line we remove the divergent part l/(2E q ) coming from vacuum. The boson energy is E h q = q 2 + m 2 . 
As lowest order contribution, the role of the DS equation for fermions is approximately equivalent to adding a 
pseudogap term A 2 g to the condensate square A 2 in the dispersion relation of quasi-particles. A few remarks about 
the pseudogap are in order. Here we define the pseudogap as the correction to fermion selfenergy from Fig. [2Kb) at 
the static limit. Its effect looks like that the real gap square gets a correction, but it is not a real gap. Actually it has 
analytical structure from which one can compute the density of states from the difermion fluctuation (H5 , [BtJ • 

The renormalized boson mass can be determined by solving the pole equation of the full propagator at static limit 
with zero gap and pseudogap, 

detp-^cp)!^^^ = det [D-^po,?)]^ + n(p 0l p)| p=A=Aj>s=0 ] - 0, (18) 

where po is a real number with analytic extension. Note that the off-diagonal parts of the boson self-energy are 

vanishing, Uir = Uri = 0. The analytical expressions for Hrr and Tin are given in App. When setting A = 0, 
we have 

IIo(po) = ^-RB.{P0,P)\ p=A =0 = U Il(P0,P)\ p=A =0 ■ ( 19 ) 

Eq. (fl8|) is then simplified as 

[(po + Pbf - m 2 b - IIoOo)] [(po - Mb) 2 -m? b - n ( Po )] = 0. (20) 
Let us focus on the following equation 

(P0 + Mfc) 2 -mg-n (p ) = 0. (21) 

In general Eq. (|2~Tj) can be understood as a complex equation by analytic extension po Po — *2 with the width 77, 
then we have 

2 

{po + Mb) 2 - -r - m l ~ Ren (po ~ ijr) = °) 

-(po + A*6)^-In>no(po-*|) = 0. (22) 
The renormalized mass for bosons can be defined by 

m\ r = (po + Mb) 2 = m b + RcTIo(po - i-j) (23) 
with a small positive 77. If the boson is stable, then 77 = and the above becomes 

m t = (Po + Vb) 2 = ml + U (p Q ). (24) 

By setting 77 = 0, we can determine the dissociation temperature T* for bosons from the first line of Eq. f22|) with 
Po = 2m — Hb = 2(m — p) > 0. It is necessary to solve the DS equation in full consistency to get T*, which is very 
involved and we reserve it for a future study. 



V. GAP AND DENSITY EQUATION 

In this section we discuss the gap and density equations which are derived from the thermodynamical potential fl 
by taking its derivative with respect to the gap A and chemical potential p. 
The gap equation reads, 

0ft _ ( m b - M 2 _ f d 3 k tanh[eg/(2T)] _ dT 2PI 
dA ~ I Ag 2 ^ J (2tt) 3 2e% d(A 2 ) 

Here the excitation energy for fermions is given by 



2A = 0. 



(25) 



4 = V^ + A 2 , 



(26) 
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with the fermion energy E k = V k 2 + m 2 and = Ek — e/i. The 2PI effective potential T2P1 is analytically evaluated 
in App. (jBj. The density equation is obtained by taking derivative of the thermodynamic potential with respect to 
the fermion chemical potential, 



d 3 k . ,„ b dT 2PJ 



e=± J 



(27r) s^(^-^) + ^f (27) 



We use an effective Fermi momentum pp to parametrize n with the relation n = p F /(3ir 2 ). Throughout the paper 
we fix the value of the effective momentum pp — 0.86. We can define from the total fermion number density n the 
density fraction for fermions, condensed and thermal bosons, and the 2PI part as 



Po 



2uA 2 
ng 2 



"^s£/(iaif[M4)-«-4)], 

P B = -Z2J j2^Y e fB{E b k -e H i b ), 



Pr 2 

which satisfies 



e=± 

1 dT 2PI 
n dfi 



Po + Pf + Pb + Pr 2 



We can treat ml as a crossover parameter and solve A, \i and A pg with the gap equation lj25j) . the density equation 
(|27ll and the pseudogap definition ifTTjl . 



VI. NUMERICAL RESULTS 



It is well known that fluctuations may change the superconducting phase transition to be of first-order, such 
as the intrinsic fluctuating magnetic field in normal superconductors [60], or the gauge field fluctuations in color 
superconductors [H3|, due to the fact that fluctuations bring a cubic term of the condensate to the effective potential 
making the Landau theroy of continuous phase transition invalid. In our case the fluctuations of di-fermions in 
superfluids also leads to a first-order phase transition. In grand canonical ensemble with chemical potential fixed, as 
shown in Fig. [3J, the thermodynamic potential fi(A) as a function of the condensate A, has a metastable state if the 
fluctuation contribution is taken into account. In the left penal with decreasing T, in the right penal with decreasing 
TOfc, we see a sudden jump of the minimum of 17(A) from A = to a non-vanishing value of A, signaling a phase 
transition from normal phase to a superfluid one. The discontinuity in the order parameter marks the first-order 
phase transition. 

We find that there are three extrema in Q(A) near the transition point, which means that there are at least three 
solutions to the gap equation with fixed chemical potential. So satisfying the gap equation is not sufficient for the 
ground state which corresponds to the global minimum of the thermodynamic potential. 

From the density equation l(27j) . the discontinuity in the gap leads to that of the particle number density. This 
situation cannot happen in canonical ensemble with fixed particle number. We have checked with the free energy 
F = ft + fin as a function of the condensate and found no metastable state. The free energy has only one extremum, 
which is the global minimum and changes continuously at the transition point. The non-zero solution to the gap 
equation corresponds to the superfluid phase, the ground state of the system. The Thouless criterion can be employed 
to determine the transition temperature. Although the order parameter changes smoothly, the phase transition is 
still a first-order one. We can take the formation of BEC in a free bosonic gas as an example. Below the critical 
temperature, the chemical potential equals to the ground state energy (zero in non-relativistic or the boson mass 

in relativistic system), so the entropy per particle si = — ^ 



chemical potential drops with a negative slope |^ 



(T < T c ) = 0. Above the critical temperature, the 

p 

(T > T c ) < giving a non-zero entropy per particle s 2 . So 
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Figure 3: The first order phase transition due to the di-fermion fluctuations. 
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there is a nonzero latent heat T c (s 2 — s\). The same phenomenon occurs to our current case with fixed total particle 
number. 

Now we present the numerical results by solving the gap equation (|25| , the density equation ([27)1 and the pseudogap 
equation ifTTj) with fixed total particle number. From the solutions we can also determine the transition temperature 
T c at which the condensate A turns to zero. 

The results for the gap A, the fermion chemical potential p,, the pseudogap A pg and the density p as functions of 
the crossover parameter mj, are shown in Fig. 21 The parameters are set to g — 1.8, T = 0.14 and m = 0.28. We 
see that the BEC/BCS regions correspond to the small/large end of to&. In the BEC region when m& is small, the 
particle population is dominated by the BEC part po while the fermion part pf is negligible, and there is no Fermi 
surface characterized by the small fermion chemical potential p. In the BCS region when nib is large, the system 
is mainly fermionic manifested by a clear Fermi surface and large fraction of fermion population. We see that the 
pseudogap A pg is small due to the small temperature chosen. The thermal density fractions ps and pf increases 
with the growing chemical potential. The contribution from the 2PI part pr 2 is large indicating the strong interaction 
between fermions and bosons. In the figure, there are breakpoints at nib ~ 0.82 for all curves because the system 
enters the normal state, i.e. A = 0, when mb is larger. In other words, the breakpoints at m& sa 0.82 indicates that 
the transition temperature T c equals to the temperature parameter T = 0.14. When m& < 0.82, the left area to the 
break point, we have T C >T = 0.14. 

In the left panel of Fig. \E[ we present the results for T c , p (at T = T c ) and A pg (at T = T c ) as functions of m&; in 
the right panel we give the results for density fractions pb, Pf, and pr 2 all at T = T c versus nib- The behavior of T c 
is similar to A, comparing Fig. [4] and Fig. [H Although the shape of p at T = T c is like p at T = 0.14 in Fig. 01 note 
that the temperature at each point of the curve of p at T = T c versus nib is different because T c itself is running with 
m&. Comparing to the result for A pg in Fig. |4j the result at T — T c is rather large especially in the BEC region where 
T c is large. The population fraction for thermal bosons pb at T c changes into pr 2 and pf with increasing nib. Pb 
intersects with pr 2 at m& fts 0.4 meaning a transition from a bosonic system to a strongly interacting one. It intersects 
with pf at nib ~ 0.9 meaning a transition to a fermionic interacting system. The intersection of pr 2 and pf occurs 
at nib ~ 1.5 indicating a transition from a fermionic interacting system to a fermionic system. 

The temperature behavior of all quantities in Fig. [4] and [5] are given in Fig. [6] with fixed nib = 0.56 (in BEC 
regime). It is easy to understand the curve of A versus T: above a critical temperature T c the condensate turns to 
zero. There is a break point on the curve of p at T c , below which p increases with increasing T but above which it 
decreases. The pseudogap is a rising function of T. The density fraction po is proportional to A 2 , so its shape is close 
to that of A. As T increases more and more condensed bosons turn to thermal bosons, fermions and interaction part. 
The thermal boson fraction ps always grows with T. The 2PI part pr 2 is also shown whose effect is large at high T. 
Due to more sensitive response to increasing temperature for bosons than fermions, the rising feature of above T c 
can be easily understood. The decreasing of pb above T c is due to particle number conservation. 



VII. SUMMARY AND DISCUSSIONS 



We investigate the fluctuation effects of di-fermion fields in the BCS-BEC crossover in a relativistic superfluid 
within the CJT formalism. We found in grand canonical ensemble the fluctuations lead to the normal-superfluid 
phase transition of first order with metastable states. In canonical ensemble, the phase transition is a first order one 
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Figure 4: The gap A, the fermion chemical potential p, the pseudo-gap A pg (the left panel) and density fractions for BEC po, 
thermal bosons ps, fermions pF, and 2PI part pr 2 (the right panel) as functions of mj with the fixed coupling constant g = 1.8, 
the temperature T = 0.14 and the fermion mass m = 0.28. 
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Figure 5: Left panel: the transition temperature T c , the fermion chemical potential p at T = T c , and the pseudogap A pg at 
T — T c versus mt- Right panel: density fractions for thermal bosons ps, fermions pF, and the 2PI part pr 2 all at T — T c 
versus nit. The parameters are g = 1.8 and m = 0.28. 
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too but without metastable states. The order parameter changes smoothly near the transition point which ensures 
the validity of the gap equation and the Thouless criterion. 

In this work, we use the bare boson mass squared mb as the crossover parameter. We find that the system is in 

the BCS/BEC regime for large/small nib- Comparing to the results in our previous work [49] with the crossover 

2 _ 2 

parameter x = — —4^2 where negative/positive x defines the BCS/BEC regime, we find that these two parameters 
are equivalent in describing the crossover behavior. We can verify numerically that there is a one-to-one mapping 
between nib and x. In this paper, we choose rrib as it is more simple. The gap equation is derived with fixed nib 
instead of x which is regarded as a derivative parameter. 

We see that the pseudo-gap A pg and the 2PI contribution to particle density pr 2 vanish as the temperature goes 
to zero. Therefore the results at low temperatures reproduce those in the mean field approximation. This is because 
as the temperature goes to zero the interaction between thermal bosons and fermions is switched off which results 
in vanishing fluctuation contributions. At finite temperature and in the middle of the crossover, pr 2 is dominant 
indicating the strong coupling between bosons and fermions. In the BEC or BCS end the number density is dominated 
by bosons or fermions respectively, during the crossover from BCS to BEC, the system undergoes a strongly interacting 
intermediate stage. This picture is a little different from that in the mean field approach where bosons and fermions 
are transformed to each other directly. 

Since the current model is a relativistic one, the anti-particles should appear in some circumstances. In the left 
panel of Fig. |4j we see that fj, = T = 0.14 at nib ~ 0.5, which means the thermal energy is comparable to the energy 
needed for exciting a pair of fermion and anti-fermion. At the deep end of BEC region the fraction of anti-particles is 
not negligible, which is partially responsible for the large condensate. In the left panel of Fig. \E[ at nib ~ 0.6 where 
T c = fi, the appearance of anti-particles drives the critical temperature and the pseudo-gap to be large. As expected 
the fluctuation effects grow with increasing temperatures. The pseudo-gap is a monotonous increasing function of the 
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Figure 6: Left panel: the condensate A, the fermion chemical potential p, and the pseudogap A pg versus T. Right panel: 
density fractions for condensed and thermal bosons po and Pb, fermions Pf, and the 2PI part pr 2 versus T. The parameters 
are set to g = 1.8 and m = 0.28 and nib = 0.56. 
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temperature, whose magnitude is seeable at critical temperature in the BEC regime. 

In this work, we expressed the effective potential in terms of bare propagators. The role of the DS equations are 
taken by the pseudo-gap and the renormalized boson mass. There is another approach of using the full propagators to 
construct the effective potential, where the complete form of the renormalization can be implemented. Another hard 
but interesting problem to solve the DS equations with full self-consistency. All these problems deserve a detailed 
study in the future. 
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Appendix A: BOSON PROPAGATOR AND SELF-ENERGY 

The bare propagator for bosons can be found from Eq. |(5]) , 

D : 

where 



1 



P^P" +fii-mi -2fi b ip 
detOD- 1 ) V 2 AWo + - m 2 b 



detOD- 1 ) = {(p a + M6 ) 2 - (E b p ) 2 ][( P0 - - (E b ) 2 } = [pi - ( M - E b ) 2 ][p 2 - (/x + E b p )^ 

We see that D RR (Q) = D H {Q) and D IR {Q) = -D m (Q). 
The boson self-energy matrix 



(Al) 



(A2) 



are evaluated as 



(A3) 



« -Tr[f nSf R S] = -8g 2 [Tr (75 Si 1755*22) + Tr^sS^sS^i)] > 
ps iTrff/ST/S] = -8g 2 [Tr^Si^sS^) - Tr^S^sS^i)] , 

= n Jfl = o, 



(A4) 
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where 



Po - M + e-Ep 



Tr[7 5 5n7 5 522] = -TV , , 

„ 7 g Po ~ (p ~ eE p ) 2 - A 2 

Tr [75 S2 1755*21] = Tr[7 5 S'i275S'i2] 

1 



Po + <?o + M + e'E p+q 
(Po + <Zo) 2 - (A* + e'£ p+g ) 2 - 



,Tr[A;A^J 



- - T E 

The results for n(po 3 p) are 



<Z§ - (/i - e£ 9 ) 2 - A 2 (po + go) 2 - 0* - e'£ p+g ) 2 - A 2 



iiijii(po,p) = n + iii 
n>/(po,p) = n - iii 



A 2 Tr[A" e A p 7 g ]. (A5) 



where 



n (po,p) 



j3 exp 
' ^TrlA^J 



P(-Pa + e'x e l+q - eie|) 



(2tt) 3 



-Po + eie^ +9 - eie^ 



-% 2 E 

, , ew , , e ' ,(ep'+g-eiM+e'e'i£;p + g)(e^+ei/i-eei£; 9 ) 
x/F(-eie^)/ F (eie^ +? )- * 



4e e ', e e 



IIi(po,P) 



2 A 2 E e * e i 



x/^(-eie,)/F(eie; +9 ) 



d 3 
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(2tt) 3 
1 

4e e ', e e 



Tr[A^A;' + J 



exp 



(3(-po + e[e^ +q - eiejj) 



(A6) 



where the quasi-particle energy for fermions is given by Eq. f26|) . Note that in renormalized quantities / constants and 
the effective potential T2P1 the full fermion propagators are used, which amounts to taking the pseudogap contribution 
into account, thus we have 



, A - ^) 2 + A 2 + A 2 g . 
In this and the next appendix, we always imply the above replacement. 



(A7) 



Appendix B: 2PI EFFECTIVE POTENTIAL 

In this appendix we give the 2PI effective potential Y2P1 using the boson self-energy in Eq. 1|A4|) as follows, 
r 2 p/ = ~Tr{D RR n RR ) - ^Tr(D n Iln) = -Tr(D HK n ) 



= ~ T E / (^s £ dT'D RR {T',p)eP° T ' £ drU (T,p)eP° T 



. rP 



d 3 p f 13 

drD RR (P -T,p)n (r,p) 



(27T) 

f _fp_ f _fq_ v TfA e A e ' 1 exp[/?(e3 ^ + ea£ % + C j4±l ~ ei ^ )] ~ 1 
7 (27r) 3 i (2tt) 3 , 2-, L « P+«J e3M6 + e2E 6 + eie e' eie e 

e,e ,ei,e^,e2,e3 r ^ y y 



xf B (e 3 fib + e 2 E b p )f F (-e 1 e e q )f F (e 1 e e p+q )— E • — , * , ( B l) 

Z£/ p Ze p+<2 Ze g 



where we have used 

n (r,p) - Tj2MPo:P)e- poT 
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